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Abstract 

We present a study of the decays 0*^71^ and —>■ D*^K^ , where the D*^ decays into D^ir^, 

with the reconstructed in the CP eigenstates K~ and tt^tt"*' and in the (non-CP) channels 
K~'K^ , K~iT~^TT~^iT~ , and K^tt^tt^. We use an unbinned maximum likelihood fit to measure the 
signal yields. Using a sample of about 123 million T(4S') decays into BB pairs, we measure the 
ratios of decay rates 



^^^+ ^ b\b- ^ D§p^7:-) = ^ 0.021(stat)+H°5^(syst), 



and provide the first measurements of 



Kon-CP ^ Vm- nT"'"'^ J = 0-0805 ± 0.0040(stat)+no32(syst), 

D[n ^ non-CP^ ' 



and of the CP asymmetry 

B{B-^D*%,K-)-B{B+^D*%,K\ ^ ^ ^ ^ 

A%p, = ^± -— ' = -0.02 ± 0.24 Stat ± 0.05 syst . 

B{B-^D^p K-)+B{B+^D^p ^ ^ ^ 
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1 INTRODUCTION 



The study of ^D^*^^ K^*^^ decays will play an important role in our understanding of CP 
violation, as they can be used to constrain the angle 7 = ar:g{—VudVub/^cdy*h) of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix in a theoretically clean way by exploiting the interference 
between the b — cus and b ucs decay amplitudes p. In the Standard Model, in the absence of 
mixing, R'cp±l ^non—CP — 1 + '"^ i 2r cos 5 cos 7, where 

_ m- - d:L-cp/cp±K^) 

r is the ratio of the color suppressed D*^K~^ and color allowed B~ — > D*^K~ amplitudes 

(r ~ 0.1 — 0.3), and 6 is the CP-conserving strong phase difference between those amplitudes. 
Furthermore, defining the direct CP asymmetry 

_ B{B-^D*S>±K-) - B{B+^D*^p^K+) 

~ B{B'^D*^p^K-) + B{B+^D*^p^K+) ' ^ ' 

we have: ^cp-|- = =b2rsin(5sin7/(l + it 2rcos5cos7). The unknowns 5, r, and 7 can be 
constrained from the measurements of -^non-CP' ^CP±^ ^CP±' '^^^ Belle Collaboration has 
reported RI^^_cp = 0.078 ± 0.019 ± 0.009 using 10.1 fb"^ of data HJ. 



2 THE BABAR DETECTOR AND DATASET 

We present the measurement of -^non-CP' ^CP+ ^^'^ ^CP+ performed using 113 fb~^ of data 
taken at the T{4S) resonance by the BABAR detector with the PEP-II asymmetric B factory. An 
additional 12 fb^^ of data taken at a center-of-mass (CM) energy 40 MeV below the T{4S) 
mass was used for background studies. The BABAR detector is described in detail elsewhere [3]. 
Tracking of charged particles is provided by a five-layer silicon vertex tracker (SVT) and a 40- 
layer drift chamber (DCH). Their identification exploits ionization energy loss in the DCH and 
SVT, and Cherenkov photons detected in a ring-imaging detector (DIRC). An electromagnetic 
calorimeter (EMC), comprised of 6580 thallium-doped Csl crystals, is used to identify electrons and 
photons. These systems are mounted inside a 1.5-T solenoidal superconducting magnet. Finally, 
the Instrumented Flux Return (IFR) of the magnet allows discrimination of muons from other 
particles. We use the GEANT4 Monte Carlo (MC) |5 to simulate the response of the detector, 
taking into account the varying accelerator and detector conditions. 



3 ANALYSIS METHOD 

In this analysis B D*^h^ candidates are reconstructed, where the prompt track is a kaon 
or a pion. D*^ candidates are reconstructed from D*^ —>■ D^tt^ decays and mesons from their 
decays to K~tt'^, K~Tr~^7r~^TT~ , K~'ir~^TT^, 7r~7r+, and K~K'^ . The first three modes are referred to 
as "non-CP modes", while the last two as "CP modes". Reference to the charge-conjugate decays 
is implied here and throughout the text, unless otherwise stated. 

Charged tracks used in the reconstruction of D and B meson candidates must have a distance 
of closest approach to the interaction point within 1.5 cm in the transverse plane and within 10 
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cm along the beam axis. Charged tracks from the D 7r~7r~'~ decay must also have transverse 
momenta > 0.1 GeV/c and total momenta in the CM frame > 0.25 GeV/c. Kaon and pion 
candidates from all decays must pass particle identification (PID) selection criteria based on 
a neural network algorithm, which uses measurements of dE/dx in the DCH and the SVT and 
Cherenkov photons in the DIRC. 

For the prompt track to be identified as a pion or a kaon, we require that its Cherenkov 
angle {9c) be reconstructed with at least five photons. To suppress misreconstructed tracks, while 
maintaining high efficiency, events with prompt tracks with Oq > 2 standard deviations (s.d.) away 
from the expected values for both the kaon and pion hypothesis are discarded; this selection rejects 
most protons as well. The track is also discarded if it is identified with high probability as an 
electron or a muon. 

Neutral pions are reconstructed by combining pairs of photons, with energy deposits larger 
than 30 MeV in the calorimeter that are not matched to charged tracks. The 77 invariant mass is 
required to be in the range 122-146 MeV/c^. The mass resolution for all neutral pions is typically 
6-7 MeV/c^. The minimum total laboratory energy required for the 77 combinations is set to 
200 MeV for vfO candidates from 1)° mesons. Only candidates with CM momenta in the range 
70-450 MeV (which we will call soft pions, tTs) are used to reconstruct the D*^. A fit is performed 
to constrain the 77 mass to the nominal vr" mass 5 . 

The mass resolution is 11 MeV/c^ for the — > K'tt^tt^ mode and about 7 MeV/c^ for all 
other modes. A mass-constrained fit is applied to the D candidate. The resolution of the difference 
between the masses of the D*^ and the daughter candidates (AM) is typically in the range 
0.8-1.0 MeV/c^, depending on the decay mode. A combined cut on the measured and soft 
pion invariant masses and on AM is also applied by means of a defined as: 

2 _ (moo-mgo)^ {m^^-m^J^ (AM - AM)^ 
X — 2 2 2 ' 

^m^yo ^rrins ^AM 

where the mean values {ruj^o, m-,^^, AM) and the resolutions (o"m^o, cr^^^, ctaa/) are measured in 
the data. Correlations between the observables are negligible. Events with x^ > 9 are rejected. 

B meson candidates are reconstructed by combining a D*^ candidate with a high momentum 
charged track. For the non-CP modes, the charge of the prompt track h must match that of the kaon 
from the meson decay. Two quantities are used to discriminate between signal and background: 

the beam-energy-substituted mass rriEs = \J {E*'^ /2 + pj • pB^/Ef — and the energy difference 

AE = E^ — E* /2, where the subscripts i and B refer to the initial e~^e~ system and the B candidate 
respectively, and the asterisk denotes the CM frame. 

The niEs distribution for B~ — > D*^h~ signal can be described by a Gaussian function centered 
at the B mass and does not depend on the nature of the prompt track. Its resolution (about 2.6 
MeV/c^) is dominated by the uncertainty of the beam energy and is slightly dependent on the 
decay mode. AE does depend on the mass assigned to the tracks forming the B candidate, and on 
the momentum resolution. The mass hypothesis of the prompt track used to calculate AE is 
denoted by a subscript AE^, where h = tt 01 K. AEk is described approximately by a Gaussian 
centered at zero and with resolution 17-18 MeV, whereas AE^^ is shifted negatively by about 50 
MeV. B candidates with niEs in the range 5.2-5.3 GeV/c^ and with AEk in the range (-100 to 
130) MeV are selected. 

Multiple candidates are found in about 10-12% of the selected events with two and four-body 
decays and in 17% of the events with — > K~-k^it^ decays. The best candidate in each event 
is selected based on the previously defined. 
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A large fraction of the background consists of continuum (non BB) events; a powerful set of 
selection criteria is needed to suppress it. The selection is optimized to maximize the significance 
of the results. In the CM frame, this background typically has two-jet structure, while BB events 
are isotropic. We define 6t as the angle between the thrust axes of the B candidate and of 
the remaining charged and neutral particles in the event, both evaluated in the CM frame, and 
signed so that its component along the e~ beam direction is positive. | cos 0^1 is strongly peaked 
near 1 for continuum events and is approximately uniformly distributed for BB events. For the 
non-CP modes, Icos^rl is required to be < 0.9 for the K'^tt^ mode, and < 0.85 for 

^ K- 7r"*"7r+7r and K tt'^it^ modes for which the levels of the continuum background 

are higher. For the CP modes, cos 6t is required to be in the ranges (-0.9 to 0.85) and (-0.85 to 
0.8) for the — > K^K^ and tt^tt^ modes respectively. Other mode-dependent selection 

criteria are applied: events with cosOto < —0.9 {\cos9tD\ > 0.95) for — > i^~7r+7r'^7r^ and 
'TT^TT^ {D^ — > TT TT"'") modcs are rejected, where 9tD is the angle between the direction of 
the -D" in the laboratory and opposite of the direction of the kaon (pion for the — > tt~tt~^ mode) 
from the in the rest frame. Furthermore the momentum of the B candidate in the CM frame 
is required to be > 0.22 GeV/c. Finally, to reduce combinatorial background in the — > K~tt~^tt^ 
final state, only those events that fall in the enhanced regions of the Dalitz plots, according to the 
results of the Fermilab E691 experiment [HI, are selected. The reconstruction efficiencies, based on 
MC simulation, are reported in Tabled 

According to the simulation, the main contributions to the BB background for B~ D*^h~ 
events originate from the decays B~ — > D^*^^ p~ and B^ — > D*~h~^ . 

An unbinned maximum-likelihood (ML) fit is used to extract yields from the data for six 
candidate types, signal, continuum background and BB background, for each choice of prompt 
track in the candidate decays B D*^h^ . The fit is performed independently for each decay 
mode. 

Three quantities from each selected candidate are used as input to the fit: AEk, w-eSj and 
the 6c of the prompt track. The distributions of AEk and rriES for the six candidate types are 
parameterized to build the probability density functions (PDFs) which are used in the likelihood 
fit. 

Correlations between the mss and AEk variables for signal events are about 5% according to 
the simulation. To account for these, signal MC events are used to parameterize the signal PDFs 
using a method based on Kernel Estimation [Jj , which allows the description of a two-dimensional 
PDF. The shapes of MC and data distributions of these variables are in excellent agreement, 
although the central values are slightly shifted (more clearly for the m-ES distribution) by different 
magnitudes for the different modes. Hence the rriES and AEk values of the signal MC events used 
for the PDF for the signal are shifted accordingly before fitting the data. Systematic uncertainties 
associated with the statistical errors on the shifted parameters are included in the final results. 

The TTi-Es PDFs for continuum background are obtained from off-resonance data with the stan- 
dard selection criteria applied. The mES distributions are parameterized with an ARGUS threshold 
function jSj: fimEs) ~ "i_E5\/l — exp [— ^(1 — y"^)], where y = niEs/'mo and uiq is the mean 
energy of the beams in the CM frame. The AEk PDFs for background candidates from the contin- 
uum are well parameterized with exponential functions whose parameters are determined by fitting 
the AEk distributions of the selected B~ D*^h~ sample in the off-resonance data. Both the mES 
and the AEk PDFs for the continuum background are taken to be the same for B~ D*^tt~ and 
B~ D*^K~ decays. MC and data distributions of mss and AEk obtained with looser selection 
criteria, in order to increase the statistics, agree well for B~ D*^Tr~ and B~ D*^K~ decays. 
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validating this assumption. The PDFs used for the CP modes are the same as for the —* K~tt^ 
mode, as very few events from off-resonance data pass the selection criteria for these modes. This 
assumption is validated by comparing the distributions of tues and AEk for the — > K~tt'^ and 
K~K^ or vr^vr^ modes obtained using a looser event selection. 

The correlation between m-ES and AEk for the BB background is taken into account with 
a two-dimensional PDF determined from simulated events, in a similar way to that used for the 
signal. 

The FID FDFs for the kaon and pion hypotheses of the prompt track are determined from 
distributions, in bins of momentum and polar angle, of the difference between the reconstructed 
and expected 9c of kaons and pions from decays in a control sample which exploits the decay 
chain D*^ — > D^tt^, — > K^tt'^ to kinematically identify the tracks. 

Initial PDFs are parameterized for each candidate type as described above. These do not 
describe exactly all the distributions of the observables used in the fit and all their correlations. 
This is in part caused by residual particle misidentification of the prompt tracks due to long and 
small non-gaussian tails in the 9c residual distributions not removed by the prompt track selection. 
Hence, we use corrected PDFs for each candidate type which are weighted sums of all the initial 
PDFs. The weights are determined by fitting pure samples of simulated signal events and of 
background from off-resonance real and MC data. The corrections affecting the signal yields are 
typically of order 1%. The fractional systematic uncertainties for the signal yields associated to 
these corrections are in the range 0.1-6.0% depending on the decay mode. 

The likelihood C for the selected sample is given by the product of the final PDFs for each 
individual candidate and a Poisson factor: 

p-N'(TvrnN ^ ^ M ■ 



i=i j=i 

where is the total number of events, Nj are the yields for each of the previously defined six 
candidate types, and N' = J2j=i 'Pj{''^ES, AEk, 9c) is the probability to measure the particular 
set of physical quantities {m-E'^,A.Ex,9c) in the i*^ event for a candidate of type j. The Poisson 
factor is the probability of observing total events when A^' are expected. The quantity C is 
maximized with respect to the six yields using MINUIT The fit has also been performed on 
luminosity weighted MC and high statistics toy MC events and it has been found to be unbiased. 
The yields thus found are corrected to account for small differences in resolutions for AEk and 
rriES between data and simulation in the parameterization of the signal. 

The preliminary results of the fit are reported in detail in Table ^ These yields are used to 
determine the CP asymmetry parameters. We measure: 

h0.007/ 



R*^p^ = 0.088 ± 0.021 (stat)Io;oo5(syst), 

h0.0039/ 
-0.0032 V 



Rl^^_cp = 0.0805 ±0.0040(stat)+°;°°^^(syst). 



Rhp+/Kon-CP = 1.09±0.26(stat)l°:J°(syst), 

A*cp+ = -0.02±0.24(stat) ±0.05(syst). 

Figure n shows the distributions of AEk for the combined non-CP and CP modes before and 
after enhancing the B D*^K component. This is accomplished by requiring that the prompt 
track be consistent with the kaon hypothesis and that mss > 5.27GeV/c^. The AEk projections 
of the fit results are also shown. 
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Table 1: Results of the yields from the ML fit. For the CP modes the results of the fit separately for 
the and B~ samples are also quoted. Errors are statistical only. The efficiencies (e) according 
to MC simulation are also reported. 



V mode 


N{B D*^tt) 


N{B D*^K) 


e (%) 




2502 ± 55 


218 ± 17 


17 




3105 ± 66 


231 ± 20 


6 




2984 ± 63 


230 ± 21 


10 


K-K+ 


245 ± 18 


21.4 ±5.2 


15 


K-K+ [B+] 


117 ± 13 


11.6 ±3.9 


15 


K'K+ [B ] 


126 ± 13 


9.8 ±3.7 


15 


7r~7r"'" 


115 ± 14 


7.4 ±4.6 


13 


7r-7r+ [B+] 


67± 11 


1.8 ±3.0 


13 


7r-7r+ [B \ 


46 ±9 


5.4 ±3.5 
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The ratio of the decay rates for B D*"'it^ and B~ — > D*"K~ is separately calculated for 
the different decay channels and is computed with the signal yields estimated with the ML 
fit and listed in Table ^ The resulting ratios are scaled by a correction factor of the order of a 
few percent, which is estimated with simulated data, and takes into account small differences in 
the efficiency between B~ — > D*^K~' and B~ D*^-k~ event selections. The results are listed in 
Table m 

Table 2: Measured ratios for different decay modes. The first error is statistical, the second is 
systematic. 



B- D*^h- Mode 


B{B^D*'^K)/B{B^D*'^Ti) (%) 


L>° ^ K--K+ 

L>° ^ K-TT+TT+TT- 


9.10±0.74l^;|^ 
7.44 ± 0.65l[]:| 
7.73 ± 0.7ll[5Jo 


Weighted Mean (non-CP) 


8.05 ± 0.40+^:i^^ 


L>o ^ K-K+ 
L>° ^ 7r-7r+ 


9.0 ± 2.31^:;^ 

7.5±4.8ij'^ 


Weighted Mean (CP) 


8.8 ±2.1+;;:^ 



4 SYSTEMATIC STUDIES 

The sources of systematic uncertainties for the yields have been identified and their contributions 
(for the measurement of -R*non-)Cp) reported in TableOl Uncertainties of the signal parameteri- 
zations of AEk and m-Es arise from the assumed shapes of the PDFs and discrepancies between real 
and simulated data; all the parameters of the AEx and ?ties PDFs have also been varied according 
to their statistical uncertainties (one s.d.) and the variations in the yields are taken (with their 
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Figure 1: Distributions of AE^ in the B D*^h sample, for K'lr'^, K~7r+7r°, A'~7r+7r+7r" 

(top, (a), (b)) and K~ ,tt~tt'^ (bottom, (c), (d)), before (left, (a), (c)) and after (right, 

(b), (d)) enhancing the B D*^K component by requiring that the prompt track be consistent 
with the kaon hypothesis and rriES > 5.27 GeV/ c^. The B~ — > D*^Tr~ signal contribution on the 
right of each plot is shown as a dashed line, the B~ D*^K~ signal on the left as a dotted line, 
and the background as a dashed-dotted line. The total fit with all the contributions is shown with 
a thick solid line. 



signs) as systematic uncertainties. For the BB and continuum background, the systematic uncer- 
tainties due to limited statistics of the MC and off-resonance data have been calculated varying the 
AEk and mss PDF by their statistical uncertainties. There are several contributions to the PID 
systematic uncertainty for the prompt track: the uncertainty due to limited statistics is calculated 
by varying each parameter of the PDF in each bin in momentum and polar angle by its uncertainty 
(keeping constant all other parameters in the same bin and all parameters in all the other bins) and 
summing all the contributions in quadrature; results obtained with alternative PID PDFs, which 
account for different Oc residual shapes and for discrepancies between data and simulation, are also 
included as systematic uncertainties. The systematic uncertainties due to the PDF reweighting 
procedure have been evaluated. Finally, errors associated to the efficiency correction factor are also 
included. 

Many of the systematic uncertainties for the signal yields have similar effects on the B~ 
D*^K- and B' D*°7r- events (they increase or decrease both fractions simultaneously), hence 
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their effect is reduced in deriving the systematic uncertainty for the measurement of the ratios, 
when all correlations are taken into account. Overall, the main sources of systematic uncertainties 
for the measurement of both -R(non_)(7p and ^cp+ ^'■^ characterization of the shapes 

of m-ES and AE^^-for the signal, to the characterization of the rriES PDFs for the background, 
to the particle identification, and to the uncertainty of the PDFs weighting procedure and of 
the efficiency correction factors. The systematic uncertainty for Aqp_^_ due to possible detector 
charge asymmetries is evaluated by measuring asymmetries analogous to those defined in Eq. ^ 
but for B~ D*^TT~ and —>■ D*^K~ events (the latter uniquely for the non-CP modes), 
where CP violation is expected to be negligible. Results for all modes are then combined, taking 
correlations into account. The measured asymmetry is —0.010 ± 0.012(stat)lQ QQ^(syst) and its 
maximum variation from zero up to one s.d. (0.022) is taken, conservatively, as a further symmetric 
systematic error on A'^p^. When combining the results for the different modes, all systematic and 
statistical uncertainties are considered to be uncorrelated, except for the contributions of the FID 
PDF (common to all modes) and of the detector charge asymmetry in the measurement of A^p_^, 
which are considered to be completely correlated. For the measurement of Rcp+/ ^non-CP 
systematic uncertainties have been considered to be uncorrelated; this assumption is conservative, 
and has negligible effect on the largely statistically limited final result. 



Table 3: Average systematic uncertainties for 



Systematic 


-^non-CP 


^CP 


Source 


^Kon-Cp/ -^non-CP C^o) 


AR*cp/R*cp{%) 




non-CP modes 


CP modes 


AEk (signal) 


+i.y 

-1.8 


+2.4 
-2.3 


AEKiqq) 


+0.3 
-0.4 


+1.5 
-2.2 


AEk{BB) 


+0.2 
-0.4 


+1.1 
-1.7 


TTiES (signal) 


+0.4 
-0.3 


+0.6 
-0.7 


rriEs (qq) 


+0.9 
-0.9 


+4.9 
-2.1 


rwES {BB) 


+1.5 
-1.6 


+3.4 
-3.3 


PDF Weights 


+2.7 
-2.7 


+1.0 
-1.2 


FID PDF 


+2.7 
-1.8 


+2.2 
-2.0 


£ Correction 


+1.5 
-1.5 


+2.0 

-2.0 



5 SUMMARY 

In conclusion, we have measured the ratio of the decay rates for B{B^ — > D*^K~) and B{B^ — > 
D*^i:~), with non-CP eigenstates. This constitutes the most precise measurement of this kind. 
We have also performed the first measurement of the same ratio and of the CP asymmetry ^cp+ 
for mesons decaying to CP eigenstates. These results, together with measurements exploiting 
B- D^K-, B- D^K*- and B' D*°K*- decays Olini, constitute a first step towards 
measuring the angle 7. All the results presented in this document are preliminary. Assuming 
factorization and flavor-SU(3) symmetry, theoretical calculations (in the tree-level approximation) 
predict: B{B- ^ D*^K-)/B{B~ ^ D^Ovr") ~ (K./Kd)'(/x//.)' -0.074, where Jk and/, are 
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the meson decays constants Our results accord with these predictions. 
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